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The fragment production in collisions of 48,50 Ca+ 12 C at 50 MeV/nucleon are simulated via the 
Isospin-Dependent Quantum Molecular Dynamics (IQMD) model followed by the GEMINI code. By 
changing the diffuseness parameter of neutron density distribution to obtain different neutron skin 
size, the effects of neutron skin thickness (S np ) on projectile-like fragments (PLF) are investigated. 
The sensitivity of isoscaling behavior to neutron skin size is studied, from which it is found that 
the isoscaling parameter a has a linear dependence on 5 np ■ A linear dependence between 5 np and 
the mean N/Z [N(Z) is neutron (proton) number] of PLF is obtained as well. The results show 
that thicker neutron skin will lead to smaller isoscaling parameter a and N/Z. Therefore, it may be 
probable to extract information of neutron skin thickness from isoscaling parameter a and N/Z. 

PACS numbers: 21.10.Gv, 24.10.-i, 25.70.Mn 


I. INTRODUCTION 

The neutron skin of a nucleus has been one of the 
hottest issues in nuclear physics. It is defined as the 
difference between the neutron and proton root-mean- 
square (RMS) radii: S np = (r 2 ) 1 / 2 — (r 2 ) 1 / 2 . The 
mechanism for the formation of neutron skin can be ex¬ 
pressed as follows. The proton-neutron interaction is 
stronger than the proton-proton or neutron-neutron in¬ 
teraction. Therefore, if the neutron number increases 
for neutron rich nuclei, the mean potential for proton 
becomes deeper, while the potential for neutron becomes 
shallower. Consequently, protons are more deeply bound, 
but neutrons are more loosely bound, which will form the 
neutron skin JlJ. The neutron skin thickness in neutron 
rich nuclei is crucial for studying the equation of state 
(EOS) of asymmetric nuclear matter due to its strong 
correlations with the symmetry energy, the slope and 
curvature of symmetry energy at the saturation density 
Thus neutron skin thickness and its effect in nu¬ 
clear reactions become an important research subject in 
nuclear physics. 

Using the Isospin-Dependent Quantum Molecular Dy¬ 
namics (IQMD) model, Sun et al. have investigated the 
neutron to proton ratio [R(n/p)\ of emitted fragments 
from projectile with different neutron skin thickness and 
shown that there is a strong linear correlation between 
R{n/p ) and 5 np , especially for peripheral collisions EH- 
Meanwhile, we have investigated the correlation between 
the ratio of triton to 3 He [I?(t/ 3 He)] and S np , which ex¬ 
hibits the similar linear relation [12]] • Projectile fragmen¬ 
tation is a well-established technique for rare isotope pro¬ 
ductions by many radioactive ion beam (RIB) facilities 
in the world. The projectile fragmentation is used not 
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only to study the reaction mechanism of nuclear colli¬ 
sion but also to extract some information on fundamen- 
tal physics fl3l4l5| . The production of heavy fragments 
will be affected by the neutron and/or proton density dis¬ 
tributions of the projectile nuclei. Thus studies on the 
effect of neutron skin thickness over the projectile-like 
fragment (PLF) will be very interesting. 

In this paper, the relationship between S np and aver¬ 
age N/Z, the isotope distribution and the yield ratios 
of PLF will be studied within the framework of IQMD 
model. The manuscript is organized as follows. In Sec. 
II, we briefly describe the method, i.e. IQMD model plus 
GEMINI code. In Sec. Ill, we present both the results 
and discussions. Finally, a summary is given in Sec. IV. 


II. THE FRAMEWORK DESCRIPTION 
A. Dynamical model: IQMD 

The quantum-molecular-dynamics (QMD) approach is 
a many-body theory to describe heavy ion collisions from 
intermediate to relativistic energies. A gen eral review 
about the model can be found in Refs. [Tol . E3- The 
IQMD model is based on the QMD model with explicitly 
inclusion of isospin degrees of freedom. It is well known 
that the dynamics at intermediate energies is mainly gov¬ 
erned by three parts: the mean field, two-body nucleon- 
nucleon collisions and the Pauli blocking. Therefore, it 
is necessary to include isospin degrees of freedom to the 
three parts for isospin dependent dynamics. In initial¬ 
ization of the projectile and target nuclei, neutrons and 
protons should be sampled separately in phase space, es¬ 
pecially for nuclei far away from the ^-stability line, of 
which the neutron and proton density distributions are 
much different. QMD model has been widely and suc¬ 
cessfully used in heavy ion collisions [Tl], 1 lSl [291 . Details 
about the description of the IQMD model can be found 

in Refs. [HUH- 
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FIG. 1: (Color online) The neutron and proton density pro¬ 
files calculated from IQMD initialization for 50 Ca at t = 0 
with f p = 1.0 but different f n . (a) f n = 1.0(5 np = 0.08 fm); 
(b) f n = 1.2(<5„ p = 0.21 fm); (c) /„ = 1.4(<5„ p = 0.37 fm); (d) 
f„ = 1.6 (5 np = 0.54 fm). 


In the present work, the following potential is used, 

U(p, r z ) = a(^) + /3( —) T + hl- r z )V c 
Po Po 2 

, rt P n Pv i rrYuk n\ 

“r C sym Tz U , (1) 

PO 

where po = 0.16 fm -3 is the normal nuclear matter den¬ 
sity. p, p n , and p p are the total, neutron, and pro¬ 
ton densities, respectively. r z is the 2 -th component of 
the isospin degree of freedom, which equals 1 or -1 for 
neutrons or protons, respectively. The coefficients cr, /3 
and 7 are parameters of the nuclear EOS. C sym is the 
symmetry energy strength due to the difference between 
neutron and proton asymmetry in nuclei. In this work, 
a = —356 MeV, /3 = 303 MeV and 7 = 7/6 are taken, 
which corresponds to the so-called soft EOS with an in¬ 
compressibility of K = 200 MeV and C sym = 32 MeV. 
V c is the Coulomb potential and U y uk is the Yukawa 
potential. Many theoretical studies show a strong corre¬ 
lation between symmetry energy and neutron skin size, 
while there is still large uncertainty of neutron skin size 
and symmetry energy [3(J [U]. As discussed in the 
reference [32j|, both the symmetry and neutron skin have 
effects on the particle productions. Therefore, it is inter¬ 
esting to study the effects of neutron skin and symme¬ 
try energy separately. The correlations between symme¬ 
try energy and the particle or fragment productions have 
also been investigated in many theoretical studies {32l - l34j . 
The flexibility of adjusting independently the size of neu¬ 
tron skin of colliding nuclei is useful for our analyses in 


this work. Moreover, in this work, different neutron skin 
thicknesses mainly influence the neutron density distri¬ 
bution in the surface region of nuclei, so our study is 
focused on the peripheral collisions. In the peripheral 
collisions, the effect of neutron skin is larger than the 
effect due to the symmetry energy [32]. 

In the initialization of IQMD, the density distributions 
of neutrons and protons are assumed to follow the Fermi- 
type form according to the droplet model (35!, 36j, 


Pi( r ) = 


p? 


1 + eX P(/ft-^4) ’ 


i = n,p, 


( 2 ) 


where p“ is the normalization constant which ensures 
that the integration of the density distribution equals 
to the number of neutrons (i=n) or protons (i= p); ti is 
the diffuseness parameter; C % is the half density radius of 
neutron or proton determined by the droplet model [36|. 

Ci = R t [1 - (bi/Ri) 2 } , i = n,p , (3) 


here bi = 0.413/if*, Ri is the equivalent sharp surface 
radius of neutron or proton. Ri and U are given by the 
droplet model. The factor is used to adjust the dif¬ 
fuseness of density distribution. In this work, f p = 1.0 is 
used in Eq.@ for the proton density distribution, while 
f n in Eq.© is varied from 1.0 to 1.6 for the neutron- 
rich projectile. Different values of S np can be deduced 
from Eq. (j2j) by changing f n . Using the density distri¬ 
butions given by the droplet model, initial coordinates 
of nucleons in the nucleus are sampled in terms of the 
Monte Carlo method. After initialization, the samples 
with satisfactory stability and expected neutron skin size 
will be selected as candidates for collisions, as described 
in Refs. [HQ- Fig. [T] shows the neutron and proton 
density profiles of 50 Ca calculated from IQMD initializa¬ 
tion. In the simulation, /„ is 1.0, 1.2, 1.4, 1.6 respectively. 
It can be seen that with the increase of /„, the neutron 
density distribution becomes more extended, while that 
of proton is almost same in the four cases. The related 
neutron skin thickness are 0.08 fm, 0.21 fm, 0.37 fm and 
0.54 fm, respectively. The initial neutron skin can last a 
long time comparable with the reaction time, although a 
small-amplitude oscillation is visible. This level of stabil¬ 
ity of the initial samples is good enough for the purpose 
of this study. As a consequence, we can study the neu¬ 
tron skin effect on the production of fragments with these 
samples. In this work, the fragments are constructed by 
a coalescence model, in which nucleons with relative mo¬ 
mentum smaller than Po = 300 MeV/c and relative dis¬ 
tance smaller than Ro = 3.5 fm will be identified as a 
cluster. 


B. Decay model: GEMINI 

The pre-fragments produced in IQMD are excited 
which are not comparable with the final products mea¬ 
sured experimentally, it is necessary to take into account 
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FIG. 2: Fragment mass multiplicity distribution in 

50MeV/nucleon 50 Ca+ 12 C with fn = 1.0 and 0.6 < 6/6 ma x < 
1.0. Open circles represent IQMD calculation at t = 200 fm/c 
and solid points are results after the evaporation by GEM¬ 
INI. The excitation energy per nucleon of the prefragment as 
a function of A pro j — Af rag calculated by IQMD simulations is 
shown in the inset. 
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FIG. 3: (Color online) The isotopic distributions plotted as a 
function of neutron number N for elements with 12 < Z < 20 
in 140 MeV/nucleon 48 Ca+ 9 Bc. Solid squares are the exper¬ 
imental data taken from Ref. [47]], and open circles are the 
calculated results with IQMD plus GEMINI. 


the evaporation effect to obtain realistic results. The 
GEMINI code is used to calculate de-excitation of these 
fragments [37|. GEMINI is a Monte Carlo code which 
allows not just light-particle evaporation and symmet¬ 
ric fission, but all possible binary-decay modes. It de- 
excites a source nucleus by a series of sequential binary 


decays until the excitation energy of the excited frag¬ 
ment being unable to go further decay. By using GEM¬ 
INI, the fragment productions are widely and successfully 
investigated with QMD model [38l - |40j . as well as, anti¬ 
symmetrized molecular dynamics (AMD) model [4ll 142] . 
In IQMD, we can construct the angular momentum and 
excitation energy for each fragment which are used as 
input to GEMINI. The angular momentum of each nu¬ 
cleon in the fragment is calculated according to classical 
mechanics 

L i = R i xP l , (4) 

where Ri and Pi are coordinate and momentum vector 
of the i— th nucleon of the fragment in the center of mass 
(CM) frame of the fragment. The total angular momen¬ 
tum of the fragment is the summation of Eq. a over 
all nucleons in it. The excitation energy of the primary 
fragment is calculated by following equation, 

Z?* _ ^excited ^ground /c-v 

— -^bind -^bind i l 5 ) 

where E^f^ ed is the binding energy of the excited frag¬ 
ment calculated from IQMD, and £bj P d nd the binding 
energy of the ground state taken from nuclear mass ta¬ 
ble [4a| . The inset of Fig. [2] displays E* as a function of 
Ap ro j — Af rag , where A is mass number and the index proj 
and frag refer to the projectile and fragment, respectively. 
Although the constructed excitation energy is model de¬ 
pendent, the excitation energy obtained from our IQMD 
calculation is comparable with that constructed from ex¬ 
perimental data and model simulations |Til ItJ . 

With the calculated excitation energy and angular mo¬ 
mentum, the fragments at t = 200 fm/c will be de- 
excited by using GEMINI. Since the main purpose of 
the present work is to study the effect of neutron skin 
thickness of the projectile on the production of heavy 
fragments, the calculations are focused on peripheral col¬ 
lisions. The reduced impact parameter is used to describe 
the centrality of collision which is defined as b/b max , 
with fo max being the maximum impact parameter. For 
peripheral reaction, the nucleon-nucleon interaction on 
the nuclear surface plays an important role, where nucle¬ 
ons are loosely bound and the density distributions are 
quite different for neutron and proton. Fig. [2] represents 
the comparison of mass versus multiplicity of fragments 
with and without GEMINI decay under the condition of 
0.6 < b/b max < 1.0. In the peripheral collision of IQMD, 
the main reaction mechanism of fragment production is 
abrasion and evaporation of nucleons or light clusters, 
which leads to the excess of heavy products and defi¬ 
ciency of intermediate mass fragments (IMFs). In con¬ 
trast, applying the afterburner significantly improves the 
fragments production since more complex fragments are 
emitted in the de-excitation calculation. 
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FIG. 4: (Color online) The isotopic yield per event plotted as function of the neutron number N with the charge number 
of fragment varying from 15 to 19 in 50 Ca+ 12 C and 48 Ca+ 12 C at 50 MeV/nucleon (upper row). The related isotopic yield 
ratios of the two reactions as function of N is displayed in the lower row. In the calculations, f n is varied from 1.0 to 1.6 for 
50 Ca+ 12 C, while it is 1.0 for 48 Ca+ 12 C. 


III. RESULTS AND DISCUSSIONS 

Before investigating the neutron skin effect in frag¬ 
ment production, the isotopic distributions calculated by 
IQMD plus GEMINI are compared with experimental 
data to evaluate the validity of the present method. In 
Fig. [3] the measured cross sections from Ref. H3] and the 
calculated results of IQMD plus GEMINI for fragments 
with 12 < if < 20 from 48 Ca+ 9 Be at 140 MeV/nucleon 
are plotted. The calculated yields are scaled to the ex¬ 
perimental data for each isotope separately. From this 
comparison, we can see that our calculations can repro¬ 
duce the shape of the isotopic distribution data quite 
well, which suggests that the IQMD plus GEMINI model 
is reasonable for calculating the heavy fragments. This 
is also demonstrated in Refs. [H, Meanwhile, Mock 
et al. have well reproduced this experimental data with 
heavy ion phase space exploration (HIPSE), Abrasion- 
Ablation (AA) and AMD [44 145| . 

As discussed in Ref. [48[, the neutron skin plays a 
significant role in the production of fragment, espe¬ 
cially for the production of neutron-rich nuclei. Hence, 
to reveal the correlation between neutron skin and the 
fragment production is interesting. Using the IQMD 
model, collisions of 48,50 Ca projectile on 12 C target at 
50 MeV/nucleon are simulated. In our simulations, the 
diffuseness parameter /„ = 1.0 is used for 48 Ca, while 
f n varies from 1.0 to 1.6 for 50 Ca. In the upper row 
of Fig. U the production yield per event for five isotopes 
with the charge number varying from 15 to 19 are plotted. 
The reason why we choose the fragments with the charge 
number varying from 15 to 19 is explained as follows. The 
main purpose of this work is to study the neutron skin ef¬ 
fects on fragments. In our IQMD, probability of fusion is 
small, so the fragments with charge number larger than 


projectile are neglected. The lighter fragments, which 
reach chemical balance at the end of the reaction will 
keep little information of the source nuclei, so the ef¬ 
fect of neutron skin thickness on these fragments is very 
small. While the heavy residue can keep as much infor¬ 
mation of the projectile as possible. Therefore, fragments 
with the charge number varying from 15 to 19 are cho¬ 
sen as a probe of neutron skin. Firstly, we compare the 
fragment isotopic distributions from reactions induced by 
projectiles with different neutron excess. By comparing 
the two systems 48 Ca+ 12 C and 50 Ca+ 12 C with f n = 1.0, 
it is demonstrated that there are more neutron-rich frag¬ 
ments produced with increasing the neutron excess of the 
projectile. This result is similar to that in Ref. [l3j]. Sec¬ 
ondly, we compare the reactions 50 Ca+ 12 C with different 
/„, which relates to different 5 np . With the increase of 
neutron skin thickness, yields of the isotope distributions 
will decrease in the neutron-rich side but have almost 
no change in the neutron-deficient side, which is con¬ 
sistent with the results by statistical abrasion-ablation 
model 0 . 

Meanwhile, it has been shown that the yield ratios of 
isotopic fragments from two similar reactions that differ 
only in the isospin asymmetry follow scaling laws [50l - [56| . 
The isotope yield ratios R 2 i(N, Z) = Y 2 (N , Z)/Yi(N, Z) 
measured in two different nuclear reactions, 1 and 2 , are 
found to exhibit an exponential relationship as a func¬ 
tion of the neutron number N and proton number Z [Hoi], 
which can be expressed by following formula: 

R 21 (N, Z) = = Cex P ( aN + P z )> (6) 

where a and /3 are two scaling parameters and C is an 
overall normalization constant. 

The isoscaling phenomenon is systematically investi- 
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FIG. 5: (Color online) The isoscaling parameter a as a function 
to 19 in 50 MeV/nucleon 50 Ca+ 12 C. 

gated for fragments not heavier than Z — 8. However, it 
would be interesting to investigate the behavior of heav¬ 
ier fragments up to the region of PLF, because the heavy 
residues, which is mainly produced by abrasion and evap¬ 
orating mechanism, may preserve some memory of the 
source configuration, for example, the neutron and pro¬ 
ton density distribution. The isotopic scaling of heavy 
projectile residues is observed in both experiment and 
theoretical simulation l53 - (58| . Consequently, from the 
research of isoscaling of heavy residues of projectile, some 
information about neutron skin of neutron-rich projectile 
could be extracted, since the neutron skin does effect the 
production of fragments [48[. This is our new start point. 

In order to investigate the neutron skin size effect on 
isoscaling, the isotopic yield ratios between two reactions 
50 Ca+ 12 C and 48 Ca+ 12 C are plotted as a function of the 
fragment neutron number, as shown in the lower row of 
Fig.H In the calculation, f n = 1.0 is used for 48 Ca, while 
f n = 1.0,1.2,1.4 and 1.6 are used for 50 Ca. From the fig¬ 
ure, we can see that larger neutron skin size will suppress 



5 np (fm) 


FIG. 6: (Color online) The average N/Z of fragment as a 
function of 5 np in 50 MeV/nucleon 50 Ca+ 12 C. 


of S np for fragments with the charge number varying from 15 


the production of the neutron-rich fragment. Although 
the neutron number and proton number of projectiles 
and targets are not changed in reactions 50 Ca+ 12 C and 
48 Ca+ 12 C, the neutron density of 50 Ca is different with 
different /„. Eq.© can be simplified for isotopic yield 
ratio with the same charge number, 

R 21 {N) =C'exp(aN). (7) 

In the grand-canonical approximation, the scaling pa¬ 
rameter a is equal to the difference of the chemical poten¬ 
tials for neutrons in the two systems, a = fj, n /T (50.f5l|. 
The parameter a is extracted by fitting ln[i? 2 i(-W) Z)\ as 
a function of N. a dependence on 5 np for isotopes with 
proton number varying from Z = 15 to Z = 19 is dis¬ 
played in Fig. [5] One can see that a decreases linearly 
with the increasing of neutron skin thickness in the pro¬ 
jectile. It is indicated that the scaling parameter dose not 
keep constant for the heavy residues with proton num¬ 
ber close to the projectile, which is different from that 
for light fragment. This different isoscaling behavior be¬ 
tween heavy and light fragments may result from differ¬ 
ent formation dynamics. The light fragments are mainly 
produced from multi-fragmentation, in which the system 
has reached chemical balance, which leads to the similar 
isoscaling behavior of light fragment. While the heavy 
residues close to projectile originate from nucleons abra¬ 
sion or evaporation. In this case, most of the nucleon 
in projectile may only act as an spectator, and accord¬ 
ingly, the system is not balanced. Nevertheless, the heavy 
residues will keep a memory of the projectile to some de¬ 
gree. From this result some information of neutron skin 
could be extracted by measuring the isoscaling behavior 
of the heavy residues. 

Finally, neutron skin effect on N/Z of heavy PLF is 
also investigated. The average N/Z for PLF with charge 
number from Z = 15 to Z = 19 as a function of S np is 
plotted in Fig. 0 It displays that N/Z decreases with the 
increasing of S np , and there is a good linear correlation 
between N/Z and 5 np . This relationship could be regard 
as another probe for neutron skin thickness. We also can 
see that for a certain S np , the N/Z will increase with 
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the increasing of charge number, which is similar to the 
results in Ref. [58] 

The neutron skin effect on fragment production, 
isoscaling, and N/Z of residues could be explained as fol¬ 
lowing. For different neutron skin thickness, the mainly 
difference of neutron density distributions is in the sur¬ 
face region. With increasing of 5 np , more neutrons are 
pushed to the surface of nuclei. Consequently, more neu¬ 
trons will be abraded in peripheral collision, which makes 
the residue less neutron-rich. Large S np makes the neu¬ 
tron become more loosely bound, and more neutrons will 
be evaporated. 

IV. SUMMARY 

Using the isospin dependent quantum molecular 
(IQMD) model followed by GEMINI, the neutron skin 
effect on the production of projectile-like fragments are 
investigated in peripheral collisions of 50 MeV/nucleon 
50 Ca+ 12 C and 48 Ca+ 12 C. By changing the neutron dif¬ 
fuseness parameter of 50 Ca to obtain different neutron 
skin thickness ( S np ), the dependence of isotopic distri¬ 
butions, neutron to proton ratio (N/Z), and isoscaling 
behavior of PLF on S np are studied. It is demonstrated 
that larger S np suppresses the production of neutron-rich 


PLF. This is because projectile with larger 5 np prefers to 
produce more free neutrons and neutron-rich light clus¬ 
ters. The isoscaling behavior in 50 Ca+ 12 C and 48 Ca+ 12 C 
are also investigated. The extracted isoscaling parameter 
a decreases linearly with the increase of S np . The depen¬ 
dence of neutron to proton ratio (N/Z) of PLF on 5 np 
displays the similar trend. With the increase of 8 np , more 
neutrons will distribute in the nuclear surface. Conse¬ 
quently, they will be abraded in peripheral collision more 
easily. In conclusion, the isotopic distributions, isoscal¬ 
ing parameter a , and neutron to proton ratios N/Z of 
PLF have dependence on 8 np . This dependence could 
probably be used to extract some information of neutron 
skin from experiments. 
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